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ABSTRACT 

We present a low energy expansion of the Kramers-Heisenberg formula for atomic hy- 
drogen in terms of (w/u>i), where uii and u> are the angular frequencies corresponding 
to the Lyman limit and the incident radiation, respectively. The leading term is pro- 
portional to (uj/lui) 4 , which admits a well-known classical interpretation. With higher 
order terms we achieve accuracy with errors less than 4 % of the scattering cross sec- 
tions in the region lo/lji ^ 0.6. In the neighboring region around Lya {lo/loi > 0.6), 
we also present an explicit expansion of the Kramers-Heisenberg formula in terms of 
Aw = (lu — uJLya)/^>Lya- The accuracy with errors less than 4 % can be attained for 
ijjjuji 0.6 with the expansion up to the fifth order of Aw. We expect that these for- 
mulae will be usefully applied to the radiative transfer in high neutral column density 
regions, including the Gunn-Peterson absorption troughs and Rayleigh scattering in 
the atmospheres of giants. 
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1 INTRODUCTION 

Hydrogen is the most abundant element in the universe and 
therefore one may often encounter astronomical situations 
associated with the radiative transfer in a region with a very 
high neutral hydrogen column density Nhi- The extended 
atmosphere around a giant star is such an example, where 
near UV photons can be significantly scattered by atomic 
hydrogen (e.g. Isliker, Nussbaumer & Vogel 1989). Another 
example may be found in searches for the first objects that 
are responsible for the reionization of the universe. When 
the universe is still partially neutral before the completion of 
the reionization process, radiation should go through regions 
with a very high neutral hydrogen column density Nhi- 

In particular, the interactions with low energy elec- 
tromagnetic waves with angular frequency lu are simply 
Rayleigh scattering, where the scattering cross section is 
known to be proportional to lu 4 in the limiting case where 
lu is much smaller than the angular frequency uiL ya corre- 
sponding to the Lya transitions. 

In quantum mechanics, the Rayleigh scattering process 
is described by a second-order time dependent perturbation 
theory, where the scattering atom suffers level transition 
twice, one associated with the annihilation of the incident 
photon and the other associated with the creation of the 
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scattered photon. The scattering cross section is known as 
the Kramers-Heisenberg formula, which is obtained by com- 
bining an infinite sum over all the bound np states with 
the energy E — —Eo/n 2 and an integral over all continuum 
states n'p with the energy E — Eo /n' 2 , where Eo is the Ryd- 
berg energy (e.g. Sakurai 1967). Since the wave functions are 
analytically known for a single electron atom, the Kramers- 
Heisenberg formula for hydrogen can be written explicitly 
in a closed form. 

However, the Kramers-Heisenberg formula is unwieldy 
due to the presence of the infinitely many atomic levels con- 
tributing to the cross section. In the case of a hydrogen 
atom in the ground state interacting with incident radia- 
tion with lu much less than Luiya, this inconvenience can be 
overcome by expanding the Kramers-Heisenberg formula in 
terms of u)/uJLy a - Devoid of any resonance in the red re- 
gion of Lya, the scattering cross section is a well-behaved 
monotonic function of lu, and the leading term is propor- 
tional to lu 4 , which admits an immediate classical interpre- 
tation. Because lu 4 dependence is the limiting behaviour of 
u)/L0Lyn <C 1, inclusion of higher order terms will be useful 
to obtain more accurate cross section values in the red vicin- 
ity of Lya. However, very near the Lya resonance, the cross 
section is well approximated by a Lorentzian and hence a 
polynomical approximation becomes poor. Lee (2003) intro- 
duced an expansion of the Kramers-Heisenberg formula near 
Lya, by computing the deviation from the Lorentzian. 

In this paper, we will provide and compute the accuracy 
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of the expansions that provide the Rayleigh scattering cross 
section redward of Lya. 



2 CALCULATION 

2.1 The Kramers-Heisenberg Formula 

The interaction of electromagnetic waves with an atomic 
electron is described using the Kramers-Heisenberg formula 
that is obtained from the fully quantum mechanical second- 
order time dependent theory (e.g. Sakurai 1969, Merzbacher 
1970). In terms of the matrix elements of the dipole operator, 
the Kramers-Heisenberg formula can be written as 
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where m e is the electron mass, ro = e 2 /m e c 2 is the classical 
electron radius, and e a ,e a are polarization vectors of inci- 
dent and scattered radiation respectively. Here, u>ia is the 
angular frequency between the intermediate state I and the 
ground state A = Is. The intermediate state / includes all 
the bound np and free n'p states, and therefore the summa- 
tion notation should be interpreted as a sum over all bound 
np states plus an integration over all continuum n'p states. 

The Wigner-Eckart theorem allows one to separate the 
matrix elements of the rank 1 tensor operator p into the 
angular part and the radial part, where the radial part is 
given by the reduced matrix elements < / || p \\ i >. (e.g. 
Merzbacher 1970). Applying the Wigner-Eckart theorem, we 
obtain 
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The dot product of the polarization gives rise to the same 
dipole type angular distribution as that for Thomson scat- 
tering or classical Rayleigh scattering (e.g. Lee & Ahn 1998, 
Lee & Lee 1997, Schmid 1989). 

Blueward of Lya the cross section shows singular be- 
haviours at resonances with bound excited states, which can 
be avoided introducing damping terms associated with the 
finite life times of excited states. However, redward of Lya, 
there is no resonance and the cross section is a well-behaved 
monotonic function of u>i. We have u < ujia for any inter- 
mediate state I redward of Lya, which allows the expansion 
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Substituting these relations into Eq. (1), we have 
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Here, we used a commutation relation 
< A\p\I >= im e Lu IA < A\x\I >, 



(4) 
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which results from the commutation relation p = 
m[x, Ho]/ih with Ho being the Hamiltonian. The leading 
term is proportional to co 4 , which is well-known result for 
the Rayleigh scattering cross section in the low energy limit. 

By performing the angular integration for unpolarized 
incident radiation and averaging over the polarization for 
an outgoing radiation, we have Jdfije'™ ' • e* Q ' 
define the dimensionless angular frequencies 
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where loi is the angular frequency corresponding to the 
Lyman limit. From this we have uj np i s = 1 — n~ 2 and 
&n' P is — 1 + ( n ')~ 2 f° r a bound np state and a contin- 
uum n'p state, respectively. In a similar way, we define a 
dimensionless position operator 
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where ao = Ti 2 /(rn e e 2 ) = 0.53 A is the Bohr radius for hy- 
drogen. 

With these definitions we may write the total scattering 
cross section 
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Here, a T = 87rrg/3 = 0.665 x 10 -24 cm 2 is the Thomson 
scattering cross section. We note that the leading term is 
proportional to cu 4 , which is a well-known result for Rayleigh 
scattering. We are interested in the higher order terms of 
uj/uji for better approximation redward of Lya. 

The matrix elements of the dipole operators < np || f \\ 
Is >,< n'p || x || Is > are easily found in textbooks on 
quantum mechanics (e.g. Berestetski, Lifshitz & Pitaevskii 
1971, Saslow & Mills 1969, Bethe & Salpeter 1967). The 
matrix elements are given by 
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for the bound states. For the continuum states, the corre- 
sponding values are given by 
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Here, the normalization condition for continuum wavefunc- 
tions is 

/>oo 

J R n ' lP {r)R n ' 2P r 2 dr = - n' 2 ), (11) 

where 8(n') is the Dirac delta function. 

With these matrix elements we may divide the terms 
involving bound states and continuum states and write the 
scattering cross section as a sum of two power series of u)/u)i 




Here, the coefficients a p ,b p are given by 
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We rewrite Eq. (12) as a single power series 




which forms a low energy expansion of the Kramers- 
Heisenberg formula. The coefficients c p are obtained from 
the coefficients a p and b p by the relation 

v 

Cp = + b q)( a p-q + bp-q)- (15) 

9=0 

In Table 1 , we listed the numerical values of the coefficients 
a p , b p and c p up to p = 9. It is noted that the continuum con- 
tribution represented by b p is small for high orders, but not 
negligible for low orders and must be included for accurate 
calculations. 



Figure 1. The cross section of Rayleigh scattering for incident 
radiation with angular frequency OAlui < ui < 0.6^;, where u>i is 
the angular frequency corresponding to the Lyman limit. The re- 
sult obtained from the full Kramers-Heisenberg formula is shown 
by the solid line, and the dashed line represents our low energy 
expansion given in Eq. (14) up to p = 9 (see the text). The crosses 
mark the values given in a tabluar form by Gavrila (1969), which 
excellently fall on the solid curve. The polynomial fit (Eq. 16) by 
Ferland (2001) to the work of Gavrila (1969) is depicted by the 
dotted line. The dashed line is always located under the solid line, 
and such a behaviour is absent in the case of the polynomial fit 
given by Ferland, for which it crosses the solid line at lj/lji ~ 0.54. 



3 RESULT 

3.1 Low Energy Expansion 

In Fig. 1, we show our result for the low energy expan- 
sion in Eq. (14) up to p — 9 by a long dashed line. In 
the figure, we also show the result from the fully quantum 
mechanical computation, which was numerically obtained 
by performing the summation and integration appearing in 
the Kramers-Heisenberg formula. The exact cross section of 
Rayleigh scattering by atomic hydrogen was presented by 
Gavrila (1967), who computed the Green function in mo- 
mentum space and provided the cross section in a tabular 
form. The crosses in the figure mark the result from the work 
of Gavrila (1967). It is noted that the cross marks are found 
on the solid line indicating excellent agreement. 

Because the Lya resonance is located at lu/loi = 0.75, 
the covergence is rather slow for lo/loi ~ 0.6, which requires 
a large number of terms in the expansion. This behaviour 
is natural because the scattering cross section changes too 
steeply near resonance for any polynomial approximation to 



be suitable. On the other hand, in the region with uj/uji < 0.5 
the agreement is almost perfect and the two results are not 
distinguished in the figure. 

For comparison purpose, we make a plot for the poly- 
nomial fit to the result of Gavrila (1967) obtained by Fer- 
land (2001). He incorporated it in his photoionization code 
'Cloudy,' in which the formula is quoted to be useful for ra- 
diation with A > 1410 A, corresponding to lj/oji = 0.647. He 
used the Lorentzian function instead for A < 1410 A, where 
the resonance of Lya dominates. The polynomial fit to the 
scattering cross section obtained by Ferland (2001) is 
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which shows the same uj 4 dependence in the leading term. 
In Fig. 1, we also show the values obtained from this poly- 
nomial fit by a dotted line. 

Because the coefficients a p ,b p and c p are all positive, 



4 Lee and Kim 



Table 1 . Numerical Values of the Coefficients a v , b v and c, 
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1.26537 


0.9157 


0.2092 


5 


81.1018 


13.5826 


0.0571 


1 


3.73766 


1.52456 


0.1368 


6 


161.896 


23.9 


0.04983 
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8.8127 


2.58886 


0.1015 
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319.001 


42.19 


0.0442 
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19.1515 


4.4587 


0.08061 


8 


622.229 


74.641 


0.03971 
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39.919 


7.75546 


0.06685 


9 


1203.82 


132.251 


0.03605 



any finite expansion of Eq. (14) gives smaller values than 
the true ones obtained from the full Kramers-Heisenberg 
formula. Therefore, the long dashed line is always located 
under the solid line. However, such a behaviour cannot be 
seen in the case of the polynomial fit given by Eq. (16), for 
which it crosses the solid line at lo/uji ~ 0.54. 



3.2 Near-Resonance Behavior 



Lee (2003) discussed the deviation of the scattering cross 
section from the Lorentzian near Lya resonance. The de- 
viation is attributed to contributions from infinitely many 
other states than 2p. According to him, near resonance, the 
Kramers-Heisenberg formula can be expanded in terms of 
Aco = iv — o»2i, where W21 = 0.75a;; is the angular frequency 
corresponding to the Lya transition. Explicitly, the expan- 
sion can be written as 
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The coefficients up to A5 are given by 
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where /12 = 


0.4162 is the oscillator strength for the Lya 


transition. 




With these coefficients, the Rayleigh scattering cross 


section near 


Lya can be written explicitly up to the fifth 


order 
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where we define u = {u> — 0.75)/0.75 = Aui/u)2i- 

In Fig. 2 we make a plot of this expansion up to the 
fifth order represented by a dotted line. We also compare 
this result with the results from the full Kramers-Heisenberg 
formula and approximations with lower order corrections. 
In this figure, the cross marks from the result of Gavrila 
(1969) fall on the solid line representing the full Kramers- 
Heisenberg formula as seen in Fig. 1. The agreement with 
the fifth order correction to the Lorentzian is excellent. 

It is very interesting to note that the approximation 



Figure 2. The cross section of Rayleigh scattering for incident 
radiation with angular frequency 0.6u>i < ui < 0.7ojj. Similarly 
with Fig. 1, the solid line represents the full Kramcrs-Hciscnbcrg 
formula and the crosses are the values appearing in Gavrila 
(1969). The dotted line represents the Lorentzian with the fifth 
order correction (Eq. 19), which gives excellent agreement with 
the full Kramers-Heisenberg formula. In the code 'Cloudy' devel- 
oped by Ferland (2001), a simple Lorentzian is used, which also 
is a good approximation, represented by the dot-dashed line. The 
dashed line represents the first order correction, which is poorer 
than the simple Lorentzian in this particula region. However, Lee 
(2003) showed that the first order correction is much better in 
the region nearer the Lya resonance. 



of the first order correction to the Lorentzian depicted by 
a dashed line is poorer than just the Lorentzian without 
any correction terms in the region u)/u>i ~ 0.7 represented 
by a dot-dash line. However, the first order approximation 
is excellent and better than the Lorentzian very near the 
resonance, i.e., > 0.70, which is not plotted here (see Lee 
2003). This behaviour is explained by the alternating nature 
of the series given in Eq. (19). 

As is done in the photoionization code 'Cloudy' by Fer- 
land (2001), the Rayleigh scattering cross section is approx- 
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imated by the Lorentzian near resonance uj/u>i > 0.647 and 
Eq. (16) for lo/uji < 0.647, which gives accurate results 
within errors not exceeding 5 %. In a similar way, the com- 
bination of Eq. (14) up to p — 9 for u/uii < 0.6 and Eq. (19) 
for ui/oji > 0.6 is accurate within errors less than 4 %. In 
particular, it should be emphasized that significant errors 
are found only near the boundary region dividing the two 
expansions. 



4 DISCUSSION AND OBSERVATIONAL 
RAMIFICATIONS 

In this paper, we obtained an expansion in terms of lo/uji 
of the Rayleigh scattering cross section by atomic hydrogen, 
which is applied in the low energy regime with lo/uji < 0.6. 
By combining this with another expansion of the Kramers- 
Heisenberg formula around the Lya resonance in terms of 
Auj = (lo — LULyu) /^Lya, we may have a wieldy and use- 
ful approximate formula for the Rayleigh scattering process 
redward of Lya by atomic hydrogen, which can be made ar- 
bitrarily accurate by inclusion of higher order terms directly 
calculated from the Kramers-Heisenberg formula. 

Rayleigh scattering by atomic hydrogen is important 
only in the presence of a scattering region with a very high 
neutral hydrogen column density Nhi- Such high column 
density media may be found in an extended atmosphere of 
a giant star where the mass loss process is already very im- 
portant. Isliker et al. (1989) considered the effect of Rayleigh 
scattering in binary systems containing a giant star. In these 
systems, for a given inclination and density distribution, the 
scattering optical depth is dependent on the wavelength, and 
therefore light curves differ according to the observed wave- 
length. This information may be quite important to investi- 
gate the mass loss process from a giant star. 

In their analysis, Isliker et al. (1989) presented the 
Rayleigh scattering cross section given by 
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where f\ k is the oscillator strength between Is and kp states. 
In their work, they neglected the contribution from the con- 
tinuum states. However, as is noted in the previous section, 
the contribution from the continuum states to the low en- 
ergy regime is not negligible, and hence caution should be 
exercised. 

In more than half of the symbiotic stars, Raman scat- 
tered O VI 6827, 7088 features are seen, which are formed 
via Raman scattering of O VI 1032, 1038 resonance doublet 
by atomic hydrogen. Being slightly less energetic than Ly/3, 
O VI 1032, 1038 doublet may excite a hydrogen atom that 
can subsequently de-excite to excited 2s state re-emitting an 
optical photon redward of Ha. This process was first identi- 
fied by Schmid (1989), where the scattering cross section is 
of similar order to that for Rayleigh scattering (e.g. Lee & 
Lee 1997, Nussbaumer, Schmid, & Vogel 1989). 

A very high column density media may be found in an 
early universe when the reionization of intergalactic medium 
initiated by the first objects was not completed. In this par- 
tially ionized universe, a significant extinction around Lya 
is expected, which will result in a big absorption trough 



known as the Gunn-Peterson effect (Gunn & Peterson 1965, 
Scheuer 1965). Thus far several quasars with redshift z > 6.2 
have been idenfied with Gunn-Peterson troughs (Becker et 
al. 2001, Fan et al. 2003). It appears that the H I col- 
umn density that is responsible for these Gunn-Peterson 
troughs are not sufficiently high for applications of our cur- 
rent work, but may be high enough to see the deviation from 
the Lorentzian approximation of the scattering cross section 
(see Lee 2003, Miralda-Escude 1998). Deeper IR search for 
higher redshifted objects may exhibit extremely high neu- 
tral hydrogen column density, where an accurate calculation 
of the cross section is required. 
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